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ABSTRACT. 

Aaplitude  scintillation  records  have  been  obtained  of  the  136  MHz 
beacon  signal  of  the  polar  orbiting  satellite  NIMBUS  4.  The 
records  have  been  aade  at  Narssarssuaq,  Greenland,  and  cover  the 
local  tiae  interval  22-02  of  the  auroral  oval.  The  average  be- 
haviour of  scintillation  vorsus  zenith  and  aziauth  is  coapared  to 
a theoretical  index  deduced  using  weak  scattering  theory.  The 
ionospheric  irregularities  are  represented  by  a power  law  spectrun 
«k"4,  where  k is  the  irregularity  wave  nuaber.  The  coaparison 
shows,  that  the  irregularities  are  elongated  along  the  field  lines 
with  an  elongation  factor  >2.5.  Probably  they  are  also  elongated 
in  the  aagnetic  eastwest  direction  with  an  elongation  factor 
slightly  greater  than  one.  The  night  tine  scintillation  corrected 
for  geoaetric  variation  shows  a aarked  increase  when  aoving  froa 
sub-auroral  to  auroral  oval  latitudes. 


1.  INTRODUCTION 


, '•  'f'v-iSS! 


Scintillation  due  to  electron  density  fluctuations  has  been  ob-  . 
served  on  line  of  sight  paths  in  the  ionosphere.  The  polar  orbit- 
ing satellite  NIMBUS  4 has  been  recorded  for  136  MHz  amplitude 
scintillation  at  Narssarssuaq,  Greenland.  The  geographic  coor- 
dinates of  the  station  are  61.2*N  and  45.4*W.  The  invariant 

latitude  is  69°.  The  satellite  is  at  a constant  altitude  of 
1100  ka.  These  data  are  -unique  because  aost  other  observations 

in  the  auroral  zone  of  polar  orbiting  satellites  have  been  made 

* 

at  lower  frequencies  and  the  recordings  at  40  and  54  MHz  frequently 
indicate  strongly  scattered  signals.  The  weak  scattering  theory  (BRIGGS 

and  PARKIN,  1963),  relating  the  scintillation  to  the  ionospheric 
irregularities,  to  a high  degree  can  be  applied  a large  percentage 
of  the  tiae  at  136  MHz.  At  40  and  54  MHz  strong,  multiple 
scattering  theory  should  be  used,  which  is  auch  aore  complicated . 


The  NIMBUS  4 records  have  been  scaled  manually  to  get  the 

. 

scintillation  indices.  All  the  night  tine  data  have  been  averaged 


to  a statistical  picture  showing  the  zenith-angle  and  propagation- 
angle  variation  of  the  scintillation.  This  variation  is  coapared 
to  a theoretical  formula  as  deduced  froa  the  weak  scattering 


theory  (COSTA  and  KELLI Y,  1976). 

In  the  deduction  of  this  fomula  the  ionospheric  irregularities 
have  boon  nodeled  by  a power- law  power  spectrum,  because  recent  in 
situ  measurements  (9YS0N  at  al.,  1974  J PHELPS  and  SAGALYN,  1976) 
have  shown  this  form  of  the  power  spectrum.  Similarly  ground 
based  measurements  ef  the  scintillation  power- spectrum  (RUFBNACH, 
1972  i SINGLETON,  1S74  } CRAN1.  197*)  indicate  that  the 


■to  -sawtjs  -c 
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irregularity  spectrum  follows  a power-law. 


2.  THE  NIMBUS-4  DATA 

The  NIMBUS-4  records  have  been  scaled  using  the  method  outlined 
by  WHITNEY  et  al . (1969).  Before  averaging  the  index  is  con- 
verted to  the  standard  S4  index.  This  has  been  done  using  a 
statistical  relationship  between  the  two  indices  as  discussed 
by  WHITNEY  (1974).  Only  nighttime  passes  from  the  local  time 
period  22-02  have  been  included.  The  hypothesis  is,  that  around 
midnight  the  irregularities  have  a constant  average  intensity 
over  a broad  range  of  the  auroral  oval.  Thus  the  main  variation 
of  the  resulting  scintillation  should  be  caused  by  the  geometry, 
which  is  expressed  through  the  theoretical  formula  to  be  discussed 
later.  Passes  obtained  during  the  period  May-December  1972  have 
been  used.  No  discrimination  between  different  levels  of  magnetic 
activity  has  been  made.  This  is  primarily  to  get  a large  data 
base,  but  further  has  the  advantage,  that  the  latitudinal  variation 
•f  irregularity  strength  is  smoothed.  The  indices  are  averaged 
in  boxes  of  10*  azimuth  x S*  zenith  angle  (Singleton,  1975) . The 
variation  of  tko  average  indices  is  illustrated  in  Figure  1.  Data 
are  missing  northwest  and  east  of  the  station  because  the  low  elsva 
tion  passes  have  not  boon  recorded.  Notice  that  because  of  the 


sun-synchronous  orbit  the  satellite  is  always  moving  fron  south 
to  north  during  the  night 


4)  the  ridge  of  enhanced  scintillation  slightly  south  of  the 
station  and  extended  in  the  magnetic  east-west  direction. 

These  sain  features  are  to  be  expected  from  the  geometry.  The 
next  step  is  therefore  to  make  a quantitative  comparison  with  the 
theoretical  index.  But  first  this  index  is  discussed  in  the  next 
section. 

3.  DISCUSSION  OF  THEORETICAL  SCINTILLATION  INDEX. 

RUFENACH  (197S)  has  described  the  mathematical  evaluation  of  the 

theoretical  scintillation  index.  When  representing  the  ionospheric 
I 

irregularities  by  a power-law  power  spectrum,  he  is  not  able  to 
find  a closed  form  for.  the  scintillation  index.  However  recently 

COSTA  and  KELLEY  (1976)  have  shown,  that  this  is  possible.  Com- 

I 

bining  their  equations  (4)  and  (S)  the  squared  S4  index  is  as 
follows:  S42  «#r  (reA)2  M*(L  sec  i)  J (1-I0  (^tc)  exp  (-^-x))  (1) 

re  ■ classical  electron  radius  ■ 2.8  x 10_15U 
A ■ radio  wave  length  ■ 2.2  m (136  MHz) 

AN  « variance  (m's)  of  ionospheric  irregularities. 
k0  ■ 2*/L0(m"1),  where  Lc  is  the  outer  scale  .of  the  ionospherio. 
irregularities. 

L seci  ■ slant  thickness  of  the  irregularity  layer, 
i • zenith  angle  of  radio  ray  at  intersection  with  irregularity 

layer. 

o • elongation  of  the  irregularities  along  the  magnetic  field 

linos., 

f*  ■ cos2*  ♦ a2»ia2f. 

ray  and  the  magnetic 

■» I ■e—UMMi 


indeed  proportional  to  Hktf/L0  in  agreement  with  (2). 

^ Z/2 

S4  is  proportional  to  \ ' as  verified  experimentally  by  CRANE 

(1976). 

In  the  case  the  irregularities  are  elongated  by  a factor  y in  the 
magnetic  east-west  direction,  a (l*0*)/83  should  be  replaced  by 

f(*»*)  ■ oy(y2cos2*  ♦ sin2*  ♦ cos2*  (cos2*  ♦ Y2sin2*)+  a2sin2*)/ 
(Y2cos2*  *a2sin2*  (y2  cos2*  ♦ sin2*))3^2  (3) 

where  * is  the  azimuth  of  the  radio  ray  in  a local  coordinate 
system  having  z-axis  along  the  magnetic  field  and  y-axis  in  the 
magnetic  east-west  direction. 

From  (2)  and  (3)  the  geometric  variation  of  is  given  by: 

S4  z/cosi  f(*,*)  (4) 

z/cosi  only  depends  upon  zenith-angle.  f(*,*)  depends  both  upon 
the  zenith  and  azimuth  angles  through  the  two  propagation  angles 
* and  *. 

Assuming  a fixed  ionospheric  height  of  3S0  km  and  using  the  con- 
stant height  of  1100km  of  NIMBUS  4,  z/cosi  increases  by  a factor 
3.6  from  zenith  to  the  horizon.  However  this  monotonic  increase 
of  s/cosi  is  partly  cancelled  by  the  variation  of  f(*,*).  If 
the  elongation  of  the  irregularities  is  larger  along  the  field 
lines  than  perpendicular  to  the  field  lines  a>Y,  f(*,*)  has  an 
absolute  maximum  for  *»0  and  decreases  monotonical ly  for  constant 
t,  that  is  — <0  for  *>0.  *«0  means,  that  the  radio  ray  propagates 

parallel  to  the  earth's  magnetic  field.  This  happens  when  the 
satellite  is  southeast  ef  the  station  and  13*  from  zenith  (com- 
pere figure  1) , Because  * reaches  larger  values  north  of  the 
station  than  to  the  south'  tower  scintillation  would  be  expected 
north  of  the  station.  However  at  largo  zenith  angles  z/cosi  takes 

■ w^§|  ill 
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over,  such  that  a valley  is  created  north  of  the  station  and  a 
flat  plateau  south  of  the  station. 

The  effect  of  introducing  y>l  is  to  extend  the  maximum  at  to 

■ 

a ridge,  which  is  oriented  in  the  direction  of  the  perpendicular 
irregularity  elongation. 

This  Tidge  nay  be  misinterpreted  as  a latitudinal  maximum  of  the 
irregularities  themselves.  This  is  analyzed  in  the  final  chapter 
where  the  relative  variations  of  the  experimental  scintillation  and 
expression  (4)  are  compared. 

4.  COMPARISON  WITH  DATA 

By  changing  the  two  elongation  factors  a and  y in  equation  (4)  it 
is  possible  to  model  the  relative  variation  of  experimental  scin- 
tillation (figure  1).  However,  the  variance  AN  (equation  1)  that 
is  the  irregularities  themselves,  may  show  some  latitudinal  depen- 
dence. To  distinguish  between  this  variation  and  the  geometric 

I 

variation,  the  latter  is  extracted  from  the  experimental  index 
¥ 

(using  (4)),  and  subsequently  the  corrected  index  is  plotted  versus 
invariant  latitude.  The  positions  of  the  invariant  latitude  circles 
are  shown  in  figune  1.  A fixed  height  of  350  km  has  been  assumed. 

Pigure  2 shows  the  latitude  profile  of  the  corrected  indices  after 
applying  three  seta  of  a,  y-velues.  To  illustrate  the  scatter  at 
a fixed  latitude  the  upper  and  lower  quartiles  of  the  corrected 
indices  are  drawn.  In  the  upper  plot  of  figure  2 (o  ■ y » 1)  only 
the  zenith-angle  variation  Ji/coti  has  been  extracted  (*(*.♦)  - 
constant,  when  a ■ y »1). 

The  profile  is  symmetric  around  a large  maximum  at  68*-69*  invariant 
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latitude.  This  is  not  an  expected  latitudinal  variation  of  the 
irregularity  strength  as  explained  in  section  two.  Further  the 
68*-69“  latitude  interval  includes  magnetic  zenith  (^*0°),  where 
the  propagation  angle  factor  f (equation  3)  has  its  maximum.  The 
profile  is  therefore  identified  with  the  propagation-angle  variation. 
From  equation  (3)  it  may  be  seen,  that  the  relative  variation  of  f 
versus  \ and  $ only  changes  little  when  a exceeds  2.5.  Here  y is 
assumed  to  be  close  to  unity.  Therefore  the  upper  plot  of  figure 
2 can  not  be  used  to  devtarnme  values  of  a,  which  exceed  2.5. 

There  is  a cone  around  magnetic  zenith,  where  the  relative  variation 
of  depends  upon  a.  However  the  resolution  of  the  data  in  figure 
1 is  not  good  enough  to  analyse  this  cone. 

Therefore  the  second  plot  of  figure  2 gives  a profile  identical  to 
profiles  using  larger  values  of  a.  It  means  the  data  should  be 
interpreted . as  follows,  that  the  irregularities  are  elongated  along 
the  field  lines  with  an  elongation  factor  a>2.5. 

However  even  this  plot  (a  « 2.5,  y * 1)  shows  a narrow  maximum  at 
68*-69*,  which  is  not  likely  due  to  a variation  of  AN. 

In  the  lower  plot  (a  ■ 2.5,  y ■ 1.3)  it  is  demonstrated,  that  it 
is  partly  possible  to  remove  this  narrow  maximum  by  putting  y-l.S. 
This  means,  that  the  irregularities  are  slightly  elongated  in  the 
magnetic  east-west  direction.  Using  larger  values  of  y create 
unrealistic  latitudinal  profiles  of  scintillation.  Notice,  that 


the  final  profile  (a  ■ 2.5,  y ■ 1.3)  of  which  all  the  geometric 
variation  has  bepn  extracted  shows  an  increase  of  scintillation  from 
60*  to  70*  invariant  latitude.  This  can  only  be  interpreted  as  a 
variation  of  the  irregularity  strength  (AN  in  equation  1)  and  is 
identified  as  moving  from  the  scintillation  boundary  to  the  auroral 


t 


oval  (AARONS,  1973).  For  convenience  the  relative  variation  of 
the  geometric  factor  (4)  with  a = 2.5,  y * 1.3  versus  zenith  and 
azimuth  is  shown  in  figure  3.  The  factor  has  been  normalized  to 
unity  at  zenith.  This  may  be  compared  directly  to  the  experimental 
data  of  figure  1. 

5.  CONCLUSIONS 

The  136  MHz  experimental  scintillation  data  shows  a relative 
variation  versus  zenith  and  azimuth  which  can  be  matched  by  the 
theoretical  index  (1)  or  the  approximate  expressions  (2)  - (3). 

Since  these  expressions  have  been  deduced  (COSTA  and  KELLEY,  1976) 
using  a power-law  ionospheric  irregularity  spectrum,  this  is  fur- 
ther evidence,  that  the  power-law  spectrum  is  realistic  also  in  the 
auroral  oval.  To  get  agreement  it  is  necessary  to  assume,  that 
the  irregularities  are  elongated  along  the  magnetic  field  lines  by 
a factor  greater  than  or  equal  to  2.5.  There  is  also  some  in- 
dication, that  the  irregularities  are  elongated  in  the  magnetic 
east-west  direction,  however  at  most  by  a factor  1.3. 

When  all  the  possible  geometric  variation  has  been  extracted,  the 
scintillation  in  the  local  time  interval  22-02  still  shows  an  in- 
crease from  60*  to  70°  invariant  latitude,  which  is  identified  as 
moving  from  polewards  of  the  scintillation  (or  irregularity)  boundary 
(AARONS,  1973). 
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FIGURE  CAPTIONS 
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Figure  1.  136  MHz  scintillation  indices  from  night  tine 

auroral  oval  plotted  versus  zenith-  and  azimuth- 
angle  of  the  transmitting  satellite.  Blank  areas 
»ean  no  data.  The  invariant  latitude  circles  62*- 
to  a fixed  ionospheric  height  of  3S0  ku. 

Figure  2.  Experimental  upper  and  lower  quartiles  of  S4 
scintillation  index  corrected  for  geometrical 
variation  and  plotted  versus  invariant  latitude. 
Three  different  sets  of  elongation-parameters  a 
(along  the  field  line)  and  y (in  the  magnetic 
•sit-west  direction)  are  used.  To  see  the  ability 
of  the  theoretical  geometrical  factor  to  match 
the  relative  variation  of  the  experimental  index 
« logarithmic  scale  is  used. 


Figure  3.  Relative  variation  of  the  theoretical  geometrical 
factor  versus  zenith  and  azimuth.  The  best  fit 
values  of  the  elongation  parameters  « • 2.S  and 
Y • 1*3  have  been  used. 


